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Abstract. We discuss quantum fidelity decay of classically regular dynamics, in 
particular for an important special case of a vanishing time averaged perturbation 
operator, i.e. vanishing expectation values of the perturbation in the eigenbasis of 
unperturbed dynamics. A complete semiclassical picture of this situation is derived in 
which we show that the quantum fidelity of individual coherent initial states exhibits 
three different regimes in time: (i) first it follows the corresponding classical fidelity 
up to time t\ ~ h^ 1 ^ 2 , (ii) then it freezes on a plateau of constant value, (iii) 
and after a time scale t% ~ minlK, 1 / 2 ^ -2 , hT 1 ^ 2 5^ 1 } it exhibits fast ballistic decay 
as exp(— const 5 4 t 2 /h~) where S is a strength of perturbation. All the constants are 
computed in terms of classical dynamics for sufficiently small effective value h of the 
Planck constant. A similar picture is worked out also for general initial states, and 
specifically for random initial states, where ii ~ 1, and ta ~ S^ 1 . This prolonged 
stability of quantum dynamics in the case of a vanishing time averaged perturbation 
could prove to be useful in designing quantum devices. Theoretical results are verified 
by numerical experiments on the quantized integrable top. 
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1. Introduction 



Squared modulus of the overlap between a pair of time evolving quantum states 
propagated by two slightly different Hamiltonians, known as the fidelity or the quantum 
Loschmidt echo, has recently attracted a lot of attention [IJ El EH IH EE EH E| • In addition 
to numerous numerical simulations, several theoretical results have been proposed to 
describe the fidelity decay in relation to the nature of the (corresponding classical) 
dynamics. Jalabert and Pastawski [2] have related fidelity decay for coherent initial 
states at very short times, namely below or around the Ehrenfest time oc log h, to the 
classical phase space stretching rate as characterized by the Lyapunov exponents. In 
more general situations, and in particular for longer time scales, fidelity decay has been 
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related to the integrated time correlation function of the perturbation through a kind 
of fluctuation-dissipation relationship |U E] . 

In a recent paper [Oj we have developed a general theory of fidelity decay based 
on a semiclassical treatment of this fluctuation-dissipation formula. It turns out that if 
the corresponding classical dynamics is fully chaotic then the decay of fidelity is, after a 
short ~ log h (Ehrenfest) timescale, independent of the structure of the initial state in 
accordance with the quantum ergodicity. On the contrary, if the corresponding classical 
dynamics is regular then the long time asymptotics sensitively depend on the structure 
of the initial state and range from a Gaussian fidelity decay for coherent initial states 
to a power-law fidelity decay for random initial states. For regular classical dynamics 
the theory |6] predicts faster decay of fidelity on a short time scale oc S^ 1 (5 =strength 
of the perturbation), as the time correlation function of the perturbation observable 
does not decay, compared to the chaotic classical dynamics where the decay time scale 
is longer oc 5~ 2 , and is longer the faster the decay of correlations we have. However, 
the fast — ballistic decay of fidelity in the case of regular classical dynamics described 
by the theory [S| does not happen in one special but important case, namely when the 
time average of the perturbation (i.e. the observable which perturbs the Hamiltonian) 
vanishes. Classically, this means that the perturbation does not change the frequencies 
of the invariant (KAM) tori in the leading order in 5, at least in the phase space region 
of interest. 

Such a case of regular classical dynamics with vanishing time-averaged perturbation 
is the subject of the present paper. Though this is not a generic case for a sufficiently 
large class of perturbations, it may emerge naturally if the system and the perturbation 
possess appropriate discrete or continuous symmetries. We will discuss general initial 
states, and specifically also coherent and random initial states. We find a very surprising 
result, namely that the quantum fidelity, after decaying for a short time (e.g. following 
the classical fidelity (BJIB1III] for coherent initial states), freezes on a plateau of constant 
value. This is purely a quantum effect and has no analogue in the classical fidelity. 
The relative time span of the plateau is of the order of inverse perturbation strength 
1/5 and can be made arbitrary large for small perturbations. However for long times 
after the plateau ends, the fidelity displays a ballistic decay with the characteristic time 
scale oc 1/5 2 , e.g. Gaussian for coherent initial states and power law t~ d for random 
initial states where d is the number of freedoms. This ballistic decay can be explained 
semiclassically due to perturbative changes of the frequencies of invariant tori in the 
second order in 5. For coherent initial states in one dimension we find and explain 
another quite surprising general phenomenon, which we call 'the echo resonance' where 
the fidelity displays sudden and significant revivals which can, under certain conditions, 
come back even to value 1. This happens at particular values of times, which depend 
on the derivative of the classical frequency with respect to the canonical action and do 
not depend on 5. 

Using the formalism of action-angle variables and its semiclassical quantization we 
derive explicit semiclassical formulae, in the leading order in h, for the fidelity in all the 
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regimes. Our results are demonstrated with high precision by numerical experiments 
using a regular quantum top which is perturbed by periodic kicking. The quantum 
saturation of fidelity, which is a central result of this paper, may also be of some practical 
importance as it provides a mechanism for stabilizing the regular quantum dynamics. 

In section 2 we define the basic quantities and study the general properties of 
the so-called echo operator whose expectation value gives the fidelity. We propose a 
useful asymptotic ansatz for the echo operator, which is used later in section 3, in 
combination with the semiclassical action-angle dynamics, to derive explicit general 
results on the echo operator and fidelity and to identify different regimes. In section 4 
we define a numerical model on which the results for two specific classes of initial states, 
namely coherent and random states, are later quantitatively validated in sections 5 and 
6, respectively. In section 7 we discuss the general picture and summarize the results. 



2. Quantum mechanics of the echo operator 

Let Hq and H$ = Hq + 5 ■ H' denote the unperturbed and the perturbed Hamiltonian, 
respectively. In order to cover the even more general case of periodically time-dependent 
(e.g. kicked) systems, say of period r, Hs(r' + r) = Hs(r'), we utilize our formalism 
in terms of the Floquet map Us = Texp(i Jq dr' 'Hs(t') /fi) , where T denotes left-to- 
right time ordered product. Dynamics is now generated by a discrete group Ug,t G Z, 
where an autonomous continuous time flow is approached in the limit r — > 0. It 
seems convenient to postulate slightly different but completely general form of a small 
perturbation 

U s = U exp(-iVT5/h) (1) 

generated by a hermitean operator V which in the leading order matches H', V = 
H' + 0(t8). We note that all the results in the paper can be trivially translated to the 
continuous time case by substituting tr — > T, r J2t'=o ~~ * Jo ^ T '■ 

Starting from the same initial state the fidelity or the Loschmidt echo F(t) is 
defined as the squared modulus of the overlap between Uq\i[)) and U^\ip), namely 

F(t) = \f(t)\\ f{t) = (4>\M s {t)\1>), (2) 

where fit) is called the fidelity amplitude and 

M s (t) = U?U* (3) 

is the echo operator. Equivalently, Mg(t) is the time-ordered propagator generated by 
the perturbation V t = U^VUq in the interaction picture |U IE] 

M s (t) = n exp (i^V^ = f exp ^ £ V^j . (4) 

The essential results on the behaviour of fidelity |Hj are then derived from the 
combination of perturbative and semiclassical considerations of the formula (J1J. For 
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example, the essential physics is contained in a linear response approximation which is 
obtained by expanding (jJJ) to second order in 5, 

-2^2 t-1 



F(t) =l~ Z A- E C(*', *") + •••, (5) 

ft t',t"=0 

C(t',t")= (V t ,V t „) - (V t ,) (V t „) , (6) 



where (•) := • is the expectation value in the initial state \ip). Thus stronger 
decay of correlations qualitatively enhances the stability of quantum motion [3J El . It 
is useful to rewrite the double-sum on the RHS of linear response formula © in terms 
of the uncertainty of the integrated perturbation operator 



t-i 



E t = r^V t/ , (7) 



namely 



F(t) = l-^{(Y? t )-^ t ) 2 } + 0(5 i ). (8) 

Here we take a slightly different route and apply the Baker-Campbell-Hausdorff 
(BCH) expansion e A e B = exp(A + B + (1/2) [A, B] + . . .) to the echo operator (gj) 

f T s *-i r 2j2 t-i t-i i 
M s (t)= exp i-E^'-^E £[K^] + ... 

I t'=0 t'=0t"=t' ) 

= exp\Ui; t 5+lr t 6 2 + ..)} (9) 



Ji V 2 

where [A, 5] := AB — BA, introducing another operator valued series 

r * = irE E [VfiVt"]- (io) 

a t'=0t"=t' 

Note that for systems with a well defined classical limit the operator T t corresponds 
to ^-independent classical observable as —i/ %[•,*] corresponds to the classical Poisson 
bracket. In the ergodic and mixing case, of say classically strongly chaotic dynamics, 
straightforward expansion of the exponential (gj) gives the Fermi-golden-rule [7] 
exponential decay F(t) = exp(— Ki) where the argument nt is precisely the double- 
sum of correlation function on RHS of (jSJ) for sufficiently long times. However in the 
opposite case of classically regular (integrable) dynamics, on which we focus in the 
following sections of this paper, the BCH form © will turn out to be particularly 
useful. 

Let us first generally discuss the expression Q from the point of view of exact 
unitary quantum dynamics. For a typical observable V, one can define a nontrivial 
time-average operator 

E 1 <_1 

V= lim — = lim - YV t >, (11) 
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which is by construction an invariant of motion, [Uo, V] = 0. In a generic case of non- 
degenerate spectrum of U , the time average is simply the diagonal part in the eigenbasis 
|n) of the unperturbed evolution, Uo\n) = e~ lipn \n), namely 

V = J £V nn \n)(n\, (12) 

n 

where V nm := (n|V|ra). In general we split the perturbation into a sum of diagonal and 
residual part 

V = V + V res . (13) 

We say that the observable V is residual if V — V res . This corresponds to ergodicity of 
this specific observable, namely V — 0, meaning that V has zero diagonal elements, and 
this is clearly a special (non-generic) situation. (j In this paper we discuss integrable 
dynamics and the class of residual perturbations. For non-degenerate eigenphases 
{(fin} the matrix elements of the second order term ()10|) in BCH expansion can be 
straightforwardly calculated in the leading order in t as 

= T E \Vnk? COt[±(^ - <p n )\ + ^(r 1 ), (14) 

tr n k 

(n|r t |m) T g-i|(¥ : 'n-< 1 5m)_|_ e -i(^n-^m)(|-t) 

1 = T( V nn-V mm )V nm , pj— ; h 0(t~ ), 71^171. 

tr h 2sm[^(ip m - ip n )\ 

Hence we see that, provided the perturbation is residual, V = 0, the limit of doubly- 
averaged perturbation defined as 

V = lim — = Z lim - y T [Vt>M (15) 

11 11 1 t'=0t"=t' 

exists and is diagonal in the eigenbasis of Uq\ 

_ _ T k + n 

V = yV nn \n)(n\, V nn = ~y \V nk \ 2 cot[|(^ - <p n )]. (16) 

n 11 k 

Note that V is again an invariant of motion, [L^V] = 0, and that, unlike for the time 
average V, its trace always vanishes trV = 0. 

In the generic case, V ^ is a non-trivial operator. For sufficiently small 
perturbation the second term in the exponent of RHS of Q can always be neglected, 
since its arbitrary (finite) norm grows as oc t5 2 following from |(n|r t |m)| < const x t [see 
eq. ()14j) ]. in comparison to the first term whose norm grows as oc td. For sufficiently 
long times, i.e. longer than the effective convergence time of the limit (jllj) . we can 
write Tit ~^ trV so the echo operator can be written as Mg(t) = exp(iV trS/h) from 
which useful semiclassical expressions for different types of initial states were derived 
[0], all showing fidelity decay on an effective time scale oc <5 _1 . In the specific case of 
residual perturbation, V — 0, the norm of the first term in the exponential on RHS of 

jj We can always choose the perturbation V to be traceless, since subtracting a constant V — > 
V — 1 tr V/ tr 1 only changes the phase of the amplitude / and does not affect the fidelity F = | / 1 2 . 
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(jHJ) does not grow in time, as we shall discuss in the next section, so the second term 
will dominate for sufficiently long times. 

Although residual perturbations are not generic in the entire set of physically 
admissible perturbations V, they may nevertheless be of particular interest in cases 
where one is allowed to shift the entire diagonal part of the matrix V nm to the 
unperturbed Hamiltonian matrix, which is diagonal by definition. Also, it is easy 
to imagine practically or experimentally important situations where vanishing of the 
diagonal part, V nn = 0, is required by the symmetry. For example, it is obvious that 
having a unitary symmetry operation R, B)R = 1, commuting with the unperturbed 
evolution, [R, Uq] = 0, and the perturbation V which has a negative 'parity' with respect 
to the symmetry operation, R)VR = —V, is sufficient to give V nn = 0. 

As the case of generic perturbations has been treated in detail in previous 
publications j3J|Hj, we shall from now on entirely concentrate on the residual case V = 0, 
unless explicitly stated otherwise. In this case we have thus found the following uniform 
approximation of the echo operator 



which is accurate, for sufficiently small 5, up to long times at least of the order 1/5 2 . 
This is a consequence of the fact, that for V — the third-order term in BCH expansion 
again grows only linearly in time ~ t5 3 , and that the fourth-order term cannot 
grow faster than ~ t 2 5 A . The rest of the paper will be dedicated to the semiclassical 
exploration of the formula (j!7j) . 

3. Semiclassical asymptotics 

Since the classical mechanics is assumed to be completely integrable (at least locally, by 
KAM theorem, in the phase space part of interest) we can write the classical limit v(j, 0) 
of the perturbation operator V in canonical action- angle variables {jk, 0k,k = 1 . . . d} 
in d degrees of freedom as the Fourier series in d dimensions 



We shall throughout the paper use lower/upper case letters to denote the corresponding 
classical/quantum observables. Note that the classical limit of the unperturbed 
Hamiltonian Ho can be written as a function ho(j) of the canonical actions only, yielding 
the well-known quasi-periodic solution of the Hamilton's equations 




(17) 



im-9 



(18) 



3t = 3, 

t = e + u(j)t (mod 2tt) 
with the dimensionless frequency vector 



(19) 




(20) 
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The classical limit of the time-averaged perturbation V is v = Vo(j) which is here 
assumed to vanish v (j) = 0. 

In quantum mechanics, one quantizes the action-angle variables using the famous 
EBK procedure ^3] where one defines the action (momentum) operators J and angle 
operators exp(im ■ 0) satisfying the canonical commutation relations 

[Jk, exp(im • 0)] = hmk exp(im • 0), k = l,...,d. (21) 

As the action operators are mutually commuting they have a common eigenbasis \n) 
labeled by <i-tuple of quantum numbers n = (ni, . . . , rid), 

J\n) = h{n + a)\n) (22) 

where < ct^ < 1 are the Maslov indices which are irrelevant for the leading order 
semiclassical approximation employed in this paper. It follows from eq. (|2~Tj) that the 
angle operators act as shifts 

exp(im • 0)|n) — \n + m). (23) 

The Heisenberg equations of motion can be trivially solved, disregarding the operator 
ordering problem in the leading semiclassical order, 

J t = jHoTt/hj e -iH Tt/h = J 

^im &t _ e iHort/h gim-0 ^-iHort/h ^ gim-u( J)t ^\m-& (24) 

in terms of the frequency operator <*>(«/). In the whole paper we use the symbol = for 
'semiclassically equal', i. e. asymptotically equal in the leading order in h. Similarly, time 
evolution of the perturbation observable is obtained in the leading order by substitution 
of classical with quantal action-angle variables in the expression (fTHj) 

V t = e ^ort/h Ve -iH rt/n ^ v m (J)e imMJ)t e im - & . (25) 

Now we are ready to write the semiclassical expressions of the two-term BCH 
expansion (|T7Jl . The operator Ej ((Zj) giving the first order BCH term can be computed 
clS cL trivial geometric series 

E t ^J2 v m (J) (l - e im ^ J >') e [m& (26) 

yielding a quasi-periodic and hence bounded temporal behaviour (due to v (j) = 0). 
Here we have introduced modified Fourier coefficients v m (j) of the perturbation 

~ , .x _ TV m (j) _ i mMj ) /2 TVjJj) (?7 , 

m{3) 2sm(§m-u>(j))- 1 ] 

As for the operator V, or T t , giving the second order BCH term, the calculation is more 
tedious. First, we plug the semiclassical dynamics (|25|) into the definition (|T5j) . Second, 
we compute the resulting commutators of the form 

[^(J)e im '■" (J)^ 'e w • e, ,l; m »(^)e im ''■ a ' (J)^ ''e im ''• e, ], 
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ih (d e 



in the leading order in K, by means of the Poisson brackets 

and substitution of variables j,e im ' by operators J,e im& . Third, we drop the terms 
for which m' + m" ^ 0, since these contain the shift operator e 1 ( T7^ '+ T7^ ")■ (^ [ e q_ (|23p]. 
giving off-diagonal matrix elements only which we know should give vanishing overall 
contribution for a residual observable, see eqs. (|14ll6|) f. As a result we find the following 
semiclassical expression 

V^v{J) % 

= -\Y.™- d i {Mj')| 2 cot (|m • . (28) 

We have derived the semiclassical expression for both terms occurring in the echo 
operator (fTTj) . for the integrated perturbation S 4 (|2l)j) . and for the doubly-averaged 
perturbation V ()28|) . However, we note that both semiclassical expressions ()26I28|) 
are subject to a potential 'small denominator' problem which is closely related to the 
one in KAM theory. This well known problem of divergence of sums over the Fourier 
index m can be avoided in a generic case. First, strict singularities at resonances 
m u = (mod 2n), where frequencies u) are evaluated in the eigenstates |n), happen 
with probability zero. Second, the near resonances give a finite total contribution if 
one assumes that the classical limit of the perturbation v(j,6) is sufficiently smooth, 
e.g. analytic in angles 6 such that the Fourier coefficients v m fall off exponentially in 
\m\. This will be assumed throughout the rest of this paper, whereas the cases of more 
singular perturbations call for further investigations. The problem is even less severe if 
the Fourier series (jl8)l is finite as will be the case in our numerical example. The results 
(|26I28|) enable us to proceed with the actual calculation of fidelity decay for 'long' and 
'short' times, where the operators V and E t are dominanting, respectively. 

3.1. Asymptotic regime of long times 

Thus for sufficiently long times t, say longer than a certain t 2 oc 5" 1 such that the 
second term Vtr5 2 /2 in BCH expansion (j9|17j) dominates the first one £ t <5, the fidelity 
(amplitude) can be written as 



f(t) = ( exp i— Vt )), for t> t 2 . (29) 




Since both operators, £j and V, have well defined classical limits, it is clear that t 2 will 
generally not depend on H, however it may depend on the precise structure of the initial 
state. Roughly it can be estimated semiclassically as 

f Actually, in the off-diagonal matrix elements of the Poisson bracket we have oscillating functions of 
time e( m +m which, for times longer than ~ l/7i, can no longer reproduce the matrix element of 
the commutator. 
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where subscript ef means effective value in the action space region of interest, i.e. where 
the initial state is distributed. Note that the actual timescale t 2 of the dominance of 
the second order can be in fact up to a factor ~ h' 1 ^ 2 longer for coherent initial states 
and for sufficiently small perturbation, as explained in section |S| 

The formula ()29)1 can be transformed, following Ref. JB], into a very useful expression 
for the semiclassical analysis by, first, writing out the average as the trace in EBK basis 
|n), second, using the fact that V is diagonal in \n) with eigenvalues = v(hn), and 
third, semiclassically approximating the sum J2 n by an integral over the action space 
h- d fd d j: 

fit) - h - d [d d jexp(i^-v(j)t)D4j), t 2 <t<t* (31) 



D^hn) : = |(^|n)| 2 . (32) 

The last step is justified for (classically long) times up to t*, such that the variation of 
the exponential in (|31|) across one Planck cell of diameter h is small, 

t* = ^V4^ ~ n ° 6 ~ 2 ( 33 ) 

\OjV\et tS 

We note a strong formal similarity between the action-space-integral (ASI) 
representation of fidelity for a residual perturbing observable (}3T| and the ASI 
representation for a generic observable jB] for times up to ~ U°5~ l 



/genericOO = h^J d d j eX P (^(j)*] A*C£ 



(34) 



This means that only vS has to be replaced by v5 2 /2 in the semiclassical analysis of 
formula (pMjl elaborated in [Oj. This shall be discussed in detail for the specific cases of 
coherent and random initial states in sections 03 and |B1 

3.2. The plateau: linear response and beyond 

For times t smaller than t 2 (|3U|) the first term in the exponential of (|17|) dominates over 
the second one, so we may write the fidelity amplitude generally as 

/(*) = (exp ( i^Et J ), for t < t 2 . (35) 



Let us first discuss the regime of sufficiently small perturbation such that the fidelity is 
close to 1, i.e. the norm of the exponential is small, ||<5£ 4 //i|| 1, so we can use the 
second order expansion of (J35|) which is precisely the linear response formula (jSJ). We 
have to compute the uncertainty of the time integrated perturbation operator E t . From 
semiclassical expression for E 4 (|2^|) we can directly compute the expectation value 

= EE ^(hn) (l - e^"^)') ^ ro V;, (36) 

where ip n := (n\ip). Similarly, we compute the expectation value of its square 

! S t) = E E v m (hn)v m ,(hn) (37) 

1 _ e ^m-u}(hn)t _ e im'-w(hn)t _|_ e i(m+m')-u}(hn)t^ ^ m 1p* + 



Quantum freeze of fidelity decay 



10 



For sufficiently many non-vanishing components ip n and/or for sufficiently large times 
t > t±, and away from certain resonance condition (all three conditions will be discussed 
in detail later) the terms with explicitly time-dependent oscillating factors exp(im • 
uj(Tin)t) can be argued to give vanishing or (semiclassically) negligible contributions. 
In fact, the time scale t\ will be determined by the condition that at typical later 
times random phase approximation can be used in dealing with the exponentials in eqs. 

Thus the above expectation values should be time independent and equal to 



their time averages 



(E t > v m ^ n _ m r n , (38) 



£ 2 ) = (£ 2 > = E | E vUKn)^~™ + E MMV'n-J ■ (39) 

n [ m^O m^O J 

This gives us a prediction that, after following a classical decay up to the short time ti, 
the fidelity should reach a constant value - a plateau and stay there up to time ti. The 
linear response value of the fidelity at the plateau is 

Further, we can easily go beyond the linear response approximation by expanding 
the formula (J33j) to all orders in 5. For this, we have to calculate the powers of the 
operator E t 

k 

E t- E l[v mi (J)e imi @ {l - e imi U:{J)t ) . (41) 

mi,...,mj.^0 1=1 

We shall use two facts in order to carry through the calculation: (i) That for leading 
semiclassical order the operator ordering is not important. (ii) Oscillatory time 
dependent terms, which in any case average to zero, typically give exponentially damped 
or semiclassically small overall contributions when used in expectation values. Thus we 
shall approximate by its time- average E fc , and the latter is calculated by selecting 
from the product (j41j) only the combinations of multi-indices which sum up to zero 
mi + • • • + mi = 0. Let us write the time average of the operator E t in terms of an 
explicit function of canonical operators 

E = 5(J,0), 

v(j,0):= Y.v m (j)e ime . (42) 
Then we calculate 

&^E k + J2{-l) l (fjE k - 1 £ S mi+ ... +rni v mi {J)---v mi (J) 

1=2 mi,...,mi^0 

= (^pE(-iy(^(^©)^/ d d x[v(J,x)]K (43) 

Please observe that is an angle-operator (which always stands in the exponential, 
so it is well defined), and x is a d— dimensional integration variable, and that in 
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order to write (J43)) we have used an integral representation of the Kronecker symbol: 
5 m = (27i)~ d f d d xe im x . Now it is straightforward to compute the power series 
SfcloO^ /fr) k '£ k /k\ by changing the summation variables to k—l and k, yielding a product 
of two exponentials 




and the plateau of fidelity (amplitude) is the expectation value of this operator 

/plateau = J ^exp \ t^(J, 0)^J exp (-^v(J , x)\ \ , (45) 

fit) = /pi a teau for t\ < t < t 2 . The very existence of such a plateau of fidelity (high 
fidelity for small 5/h) is very interesting and distinct property of quantum dynamics. 
Note that the time scale t\ only depends on the unperturbed dynamics, namely on the 
property of the operator £ 4 , so it cannot depend on the strength of the perturbation, 
ti = O(5 ). Thus the range of the plateau, i.e. t%jt\ oc can become arbitrarily large 
for small 5. The formula (J45|) becomes very useful whenever one is able to semiclassically 
compute the quantum expectation value in terms of classical phase space integrals. We 
shall present this derivation for two extremal cases of coherent and random initial states, 
in sections El and respectively. 

Note that the formulae f)44|45j) may be very useful in a general case whenever one 
has to calculate an expectation value of the form ^exp ( i-rS()> where is a time- 
integrated quasi-periodic process with a zero time-average. 



4. Numerical example: integrable top 

For numerical illustration of the above theory we take a spin system with the following 
one-time-step unitary propagator 




with parameters a and (3. Sk, k = x, y, z are standard quantum angular momentum 
operators with a fixed magnitude S of angular momentum and with the SU(2) 
commutator [S k , Si] = ie k i m S m . 

The semiclassical limit is obtained by letting S — » oo while the classical angular 
momentum hS = 1 is kept fixed, so that the effective Planck constant is given by 
h = 1/5. The classical map corresponding to the one-time-step propagator U can be 
obtained from the Heisenberg equations of angular momentum operators in the S — > oo 
limit. Defining by (x, y, z) = (S x , S y , S z ) / S a point on a unit sphere, we obtain a classical 
area preserving map: 

x t +i = x t cos (a(zt - 0j)-yt sin (a(z t - 0j) 

y t+ i = y t cos (a(z t - (3)) + x t sin (a(z t - (3)) (47) 
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This classical map represents a twist around z-axis. We note that it corresponds to the 
stroboscopic map (j!9|) with an arbitrary unit of time, so we put r = 1, of an integrable 
system with the Hamiltonian ho(j) = \ot{j — (3) 2 generating a frequency field 

^J) = ^ = a(j-f3). (48) 

Here we used a canonical transformation from a unit-sphere to an action-angle pair 
(j,9) G [-1,1] x [0,2tt), namely 



x = y 1 — j 2 cos#, y = y 1 — j 2 sin 8, z = j. (49) 

Now we perturb the Hamiltonian by periodic kicking with a transverse pulsed magnetic 
field in x direction, 

1 / 00 

h s (j, 6, t) = -a(j - (3) 2 + Syfl - j 2 cos 9 £ 5{r-k). (50) 

k=— oo 

Perturbed quantum evolution is given by a product of two unitary propagators 

U s = U exp(-i8S x ), (51) 
so the perturbation generator is 

V = SJS. (52) 
The classical perturbation has only one Fourier component, namely 

v{j,6) = y/l-j*co80, v ±1 (j) = yi-f, (53) 

whereas vq = indicating that the time- average vanishes v = 0, and V = 0. 

In our numerical illustrations two different types of initial states will be used. SU(2) 
coherent wavepackets are used to probe the correspondence with the classical fidelity, 
while random states are used to investigate the other end - states without a classical 
correspondence. The parameter a in Uq ()46|) will be always set to a = 1.1, while (3 = 
for coherent initial states, and (3 = and (3 = 1.4 for random initial states. The reason 
for choosing nonzero shift (3 for random states will be explained later. We should stress 
that we have performed calculations also for other choices of regular Uq, also in KAM 
regime, e.g. for precisely the same model and parameter values as used in Ref. jlOj . 
and obtained qualitatively the same results as for the presented case of unperturbed 
dynamics. The coherent state written in the canonical eigenbasis |m) of the operator 5* z 
and centered at the position n = (sin-$* cos<^*, sint?* sin ip*, cost?*) on a unit sphere is 

S , 9 ci xl/2 

\0*,(p*)= £ [ s ^ m ) cos(r/2) 5+m sin(^72) 5 - m e- im ^|m)- (54) 
The corresponding classical density reads ^T] 

p cl (tf, <p) = exp {-S ((* - r ) 2 + {<p- V *) 2 sm 2 <?) }. (55) 

In the numerical experiments reported below the coherent initial state will always be 
positioned at the point (#*,</?*) = (1, 1). 
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5. Semiclassical asymptotics: coherent initial state 

Let us now study an important specific case of a (generally squeezed) coherent initial 
state which can be written in the EBK basis as a general Gaussian centered around a 
phase space point (j*, 6*) 

(n\j*,0*) — ( - J |detA| 1/4 exp|- — (hn-f)-A(hn-f) -in-0*J, (56) 

with A being a positive symmetric d x d matrix of squeezing parameters. Note that the 
shape of the coherent state is generally only asymptotically Gaussian, as Ti — > 0, due 
to cyclic and discrete nature of coordinates 6 and j, respectively. Let us also write the 
structure function (I32|) of our coherent state (foT)j) 

^(i)=(-J |detA| 1/2 exp[--(j-j*)-A(j-j*)| (57) 

which is normalized as h~ d f d d jD^p(j) = 1. 

For example, for SU(2) coherent state of a quantum top (p>fj) when written in the 
asymptotic form (|56|). the squeezing parameter reads A = 1/ sm 2 {}*. 



5.1. The plateau: linear response and beyond 

In the regime of linear response, valid for sufficiently small 5, we simply evaluate 
the general expressions (|36I37|) for the particular case of a coherent initial state (|56|) . 
namely we write ip n = (n\j*,0*). First we will show that the time dependent terms 
in expectation values of the powers of S t indeed vanish for t > t± as stated in section 
13.21 We recall the assumption that the perturbation v(j,6) is sufficiently smooth, 
e.g. analytic in 6, so that the Fourier coefficients v m (j) decrease sufficiently fast, e.g. 
exponentially, or only a finite number of v m (j) is non-vanishing. What we actually need 
here is that an effective number of Fourier components is smaller than the width of a 
wavepacket which is ~ hT 1 ^ 2 . This means that within the range of Fourier series over 
m, or ml we can approximate 

C-™VWm' = D^hn)e-^ m+m '> e \ (58) 



Let us estimate the general time dependent term of expressions ()36|37|) . where all factors 
with non-singular classical limit are combined together and denoted as g(j), by means of 
expanding the frequency around the center of the packet u>(j* + x) — + + . . ., 

where Q is a matrix Qki = du>k(j*)/dji, followed by d- dimensional Gaussian integration 



n J 

fh\ d/2 r / 1 

= g(j*)e imMr)t I -J |detA| 1/2 y d d x exp f --x ■ Ax + itm ■ fix 



g(j*)e im u}{r)t exp ( -^-m • ttA- l tt T m j . (59) 
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Figure 1. The short time fidelity for a quantized top a = 1.1, /3 = is shown for 
the coherent initial state, for S — 200 (a), and S = 1600 (b), with a fixed product 
SS = 0.32 being well described by the linear response. In (c) we show S = 1600 
and stronger perturbation with 5S — 3.2. Note that the time axis is rescaled as t/t\. 
Symbols connected with dashed lines denote the corresponding classical fidelity. The 
horizontal chain line denotes the theoretical value of the plateau 1)68(1 , while the vertical 
chain line denotes the estimated theoretical value for t-i l(72|l . In (b,c) we also indicate 
fractional 2irk/p resonances with k/p marked on the figure (see text for details). 



We see that all these terms decay to zero with Gaussian envelopes with the longest time 
scale estimated as 

h= i-mm{rn-QA' 1 Q T rn)j oc h~ 1/2 . (60) 
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Note that the decay of (|59|) is absent if Q = 0, e.g. in the case of d- dimensional harmonic 
oscillator. There may also be a general problem with the formal existence of the scale t\ 
(JBOJ) if the derivative matrix Q is singular, but this may not actually affect the fidelity 
for sufficiently fast converging or finite Fourier series (|18|). 

Thus we have shown that for t > ti expectation values (j36pr>7j) are indeed given by 
time-independent expressions ()38I39|) which in the case of a coherent initial state (Jo7)j) 
evaluate to 

(s t ) = £ v m (ry m - d * =v(f,o*), (6i) 

= (v(j*,0*)f+ E I^Cf)| 2 . (62) 

The variance v = (E 2 ) — (S t ) 2 which determines the plateau in the fidelity (|SJ) is the 
second term on RHS of eq. (jo^j) . In terms of original Fourier coefficients ()27|) the final 
linear response result reads 

5 2 

F coh {t) = 1- ^coh + 0(5 4 ), t>t h 

^ T 2 \v m (j*)\ 2 , . 

^ 4sm 2 (im-a;(j*)) 

Beyond the linear response approximation, the value of the plateau can be computed 
by applying a general formula for / p i a teau (|3H|) . We shall make use of the fact that for 
coherent states we have the expectation value 

(exp(-(i8/h)g(J, 0))> S ex P (-(i<V%(j*, 0*)), (64) 

for some smooth function (7, provided that the diameter of the wavepacket ~ \ffi is 
smaller than the oscillation scale of the exponential ~ h/5, i.e. provided 5 <C h 1 ^ 2 . 
Then the squared modulus of /plateau rewrites as 



plateau — ^ 27r )2d 



J d d a;exp (— jr^U* 



,x 



(65) 



The expression for i/ coh is of course just the lowest order expansion of Fpiateau ()65|) ■ 
It is interesting to note that the angle 0* does not affect the probability -F p i ateau as it 
only rotates the phase of the amplitude /plateau- 

For smaller times t < t± <x n 1 the quantum fidelity is expected to follow the 
classical fidelity as defined by the overlap of two initially Gaussian classical phase space 
densities evolved under slightly different quasi-regular time evolutions [see Ref. |B] for 
a definition and linear response treatment of the classical fidelity]. More precisely, the 
quantum fidelity can be written in two equivalent ways as 

F(t) = (2nh) d J d d qd d pW utoM {q,p)W utw {q,p) (66) 

= (2nh) d J d d qd d pW m (q,p)W Mdim {q,p) (67) 

where W\ 7 p\(q, p) is the Wigner function of some state The corresponding classical 
fidelity is defined by the same formula if (2iih) d ^ 2 Wut\^)(q, p) is substituted by the 
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evolving classical phase space density, a solution of the corresponding classical Liouville 
equation, with the initial condition (2irh) d / 2 W\ 1 p)(q,p) which is proportional to the 
Wigner function of the initial state Of course, this only makes sense if the function 
W\^)(q,p) is strictly nonnegative so that it corresponds to some classical state, such 
as for example for a coherent state where it is a Gaussian. Indeed, time t\ ~ h^ 1 ^ 2 
may also be interpreted as the integrable Ehrenfest time up to which phase space point- 
like quantum-classical correspondence will hold. Namely it is consistent with the time 
needed for a minimal uncertainty wavepacket of diameter ~ h 1 ^ 2 to spread ballistically 
over a region of the classical size (~ h°) of an invariant torus. After this time, quantum 
wavepacket will start to coherently interfere with itself, e.g. its Wigner function will 
develop negative values, so the strict quantum-classical correspondence will stop (see 
subsection 15. 4J) . Therefore we expect initial agreement between the classical and the 
quantum fidelity up to time t\ and after that the classical fidelity of a regular dynamics 
with a residual perturbation decays with a power law oc (5h- 1/2 t)- d (factor h 1/2 comes 
from the size of the corresponding classical density) [see Refs. [HI E]] whereas the quantum 
fidelity freezes to a constant value as computed by our semiclassical theory (jBllj) . 

This picture is nicely confirmed by the numerical experiment with a quantized 
integrable top ()46|) as shown in figure ^ where we choose zero shift (3 = 0. Agreement 
between the classical and quantum fidelity up to ti ~ y/S (h = 1/S) can be nicely 
observed. After t 1 the fidelity stays constant until t 2 , the point where fidelity again 
starts to decrease. This second timescale t 2 will be discussed in the next subsection. 
The value of the plateau can be calculated specifically for our model by means of the 
semiclassical expression for -F plateau and using Fourier modes of our numerical model 
(}53|) . We get v(j, 6) = — |a/1 — j' 2 sm(9 — sin (aj/2) and the integral occurring in 

-^plateau is elementary and gives 

bateau = Jg ( 5S 2 ^^ /2) ) > ( 68 ) 

with Jo being the zero order Bessel function. Agreement with this theory is excellent 
both in the linear response regime (figure HJi,b) and also for strong perturbation (figure 
Uk). Observe also a power law decay of the classical fidelity F c \ ~ (tS/h)' 1 beyond the 
regular Ehrenfest time t > ti for strong perturbation in figure ^p. 

Actually the calculation of -F p i a teau can be generalized to any perturbation with a 
single nonzero Fourier mode ±mo with the result 

F -- T 2 ( T6 fegoCQj ^ (aq\ 

^plateau - J Q ^ ^ ^ ^ )/2} ) > ^ J > 

whereas for a more general multi-mode perturbations we have to evaluate the integral 
(|65|) numerically. 

5.2. Asymptotic regime of long times 

After a sufficiently long time t 2 of order 5 _1 the second order term in BCH expansion 
will start to dominate over the echo operator so the fidelity can be computed with our 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

x=5 2 S 1/2 t 

Figure 2. Long time ballistic decay of the fidelity for a quantized top with a = 
1.1, (3 = and coherent initial state is shown for cases S = 200 (a), and S = 1600 (b), 
of weak perturbation SS = 0.064, and for strong perturbation S = 1600, 6S — 3.2 (c). 
Chain curves indicate theoretical Gaussian (|71|l with analytically computed coefficients, 
except in the case (c), where we multiply the theoretical Gaussian decay by a prefactor 
0.088 being equal to the theoretical value of the plateau (|6^)l . and rescale the exponent 
of the Gaussian by a factor 0.8 taking into account the effect of non-small first term 
in the exponent of (|1 T|) . Note that in the limit S — > oo the agreement with the 
semiclassical theory improves and that the size of the resonant spikes is of the same 
order as the drop in the linear response plateau. The insets show the data and the 
theory in the normal scale. 

semiclassical formula ((HI]) . The straightforward calculation follows exactly the one for 
a generic perturbation in Ref. j^j, paragraph 2.2.2, where vS has to be replaced by vS 2 /2 
so we shall not repeat it here. Such a Gaussian approximation is justified provided 
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the stationary point of the exponent is not moved appreciably from the center j* of 
the packet (|57|). This implies that t5 2 <C 1 which is the same condition as required by 
substituting the sum over quantum numbers by the integral over the action space (jHTjl . 
The final result reads 

/-*(*) * exp (-( W ■ A^u) — + -^—j , « := r-^i (70) 

where the vector u is just a gradient of the classical observable v at the center of the 
wavepacket. Thus we have derived a Gaussian decay of the fidelity 

Fcoh(t) = exp j- (^-) 2 | , t coh = (u ■ A-^-^M^ (71) 

on a timescale t co h oc h 1 / 2 5~ 2 . Indeed, it can now be checked that in the semiclassical 
regime of small h the fidelity decays well before the time limit ([33)1 . t* ~ h°5~ 2 , of 
our approximations. The formula (J7T|) can only be expected to be accurate provided 
the plateau is close to 1 and hence described within the linear response approximation. 
Only in such a case can the effect of the first term £ t in the exponential of the echo 
operator (|T7j) really be neglected, in the opposite case we can correct the Gaussian (J7T|) 
by multiplying it with the plateau value and adjusting the coefficient in the exponential 
(see e.g. figure Eh)- 

— 1/2 

In such a regime of small perturbation, 5 < ^ coh n, we determine the crossover 
time ti by comparing the linear response formula ()63j) with the decay law (|71)l. namely 
1 — (^Acoh) 2 = 1 — 5 2 v C oh/fr 2 - For stronger perturbation, namely up to 5 ~ Vh, time 
scale t 2 can be simply estimated by t co h, so we have a uniform estimation 

t 2 = minil, — z^cohl^coh = niiniconst h^^d^ 1 , const Ti^^d^ 2 } (72) 
ft. 

We note that the crossover time t%, for coherent initial states and for a small perturbation 
5 < z/ co j/ 2 fr, is in fact by a factor %~ 1 ^ 2 longer than the estimate ()3()j) . This is due 
to the fact that coherent states are strongly localized in action coordinates (quantum 
numbers) for small h. Therefore, for small 5, the operator Vt5 2 , although it may already 
be dominating £ 4 <5 in norm, will only effectively rotate the overall phase of a coherent 
initial state since it is diagonal in |n) and thus will not (yet) affect the fidelity. So 
the estimate (|5U|l is expected to be valid only for initial states whose relative support 
in quantum number lattice is not shrinking as h — > 0. It is interesting that for the 

1 /2 

strongest allowed perturbation <5 max = const n 1 for our semiclassical theory to be 
valid, the estimate (fT2~J) agrees with the general estimate (JHOJ) , ti ~ l/<$- 

Timescale t 2 can be seen in figure ^ as the point of departure of fidelity from the 
plateau value. Using our model and the position of initial coherent state this can be 
calculated to be t 2 = min{0.57-\/5 /5, 0.57/(5 2 v / 5')} (similarity of numerical prefactors is 
just a coincidence) which is shown with vertical chain lines in figure ^ The theoretical 
position of t 2 is shown with a vertical chain line and is given by t co h for a strong 
perturbation 8S = 3.2 in figure Hh while it is t co hZ / coh ( V^ m figures Hk,b. The long 
time decay of fidelity is shown in figure El Theoretical Gaussian decay (JHJ), shown 
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with a chain curve, is again confirmed with an analytical formula for the decay time 
^coh = 0.57(5^ 2 ^ 1//2 evaluated at the particular position of the packet. Note that we do 
not have any fitting parameters, except in the case of a strong perturbation (5S ^> 1, 
figureEfc) where the prefactor and the exponent of a Gaussian had to be slightly adjusted 
due to the non- negligible effect of the first term in (fTTjl [see caption for details]. Quite 
prominent feature in figures and 121 are also "spikes" occurring at regular intervals, 
where the fidelity suddenly increases or wildly oscillates. These will be called the echo 
resonances and are particular to one-dimensional systems. 

We should remark that, although we obtain asymptotically Gaussian decay of 
fidelity for a single coherent initial state, one may be interested in an effective fidelity 
averaged with respect to phase space positions of initial coherent state [TO]- In such a 
case one may typically get a power law decay due to possible points in the phase space 
where the theoretical expression for t co h diverges (at the positions of zeroes of u(j*) of 
eq. (f70J)). but still on a time scale oc 5~ 2 . Note that this effective power law decay is a 
general scenario and is not particular to the case of vanishing time-average perturbation 
(in the case of V ^ the decay time scales as 

5.3. Echo resonances in one dimension 

Let us now discuss the behaviour of the fidelity for initial wavepackets in the regime of 
linear response approximation in some more detail. We shall consider possible deviations 
from the random phase approximation in the time dependent exponentials of eqs. (|36I37|) 
which have been invoked previously in order to derive time independent terms ()38J39j) 
of the fidelity plateau pPj). and also Specifically we will explain the resonances 

observed e.g. in figure [TJ 

For such a resonance to occur the phases of (|36J37|) have to build up in a constructive 
way and this is clearly impossible in a generic case, unless: (i) We have one dimension 
d = 1 so we sum up over a one-dimensional array of integers n in the action spaced . (ii) 
The wavepacket is localized over a classically small region of the action space/lattice 
such that a variation of the frequency derivative du(j)/dj over this region is sufficiently 
small. The quantitative conditions for the occurrence together with the strength and 
the shape of such resonances are discussed below. 

In this subsection we thus consider a one-dimensional case, d — 1. Again we 
study time dependent terms of ()36I37|) which can all be cast into a general form (J59*Jl . 
however, now the time is not small enough to enable the sum over quantum number to 
be estimated by an integral. On the contrary, we seek a condition that the consecutive 
phases in the exponential build up an interference pattern. 

5.3.1. 2it -resonance: Let us expand the frequency around the center of the packet 

u(j) = u,* + (j - j V + J*?"" + ■■■ (73) 

X In more than one dimension we would clearly need a strong condition on commensurability of 
frequency derivatives over the entire region of the action lattice where the initial state is supported. 
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Figure 3. Structure of echo resonances for coherent initial states of the modified 
quantum top U = exp {-\S[a{S z / S) 2 /2 + j(S z /S - j*) 3 /6]), a = 1.1, S = 200,5 = 
1.6 ■ 10 -3 , for increasing value of lo" = 7 = (a), lo" = 1 (b), and lo" — 4 (c), which 
weakens and broadens the resonances. Note that in (a), lo" = we have the same data 
as in figure^,. Vertical chain lines show theoretical times t r /2 for 7r, and t Y for 2ir 
resonance. 



where 

lo* = lo(j*), lo' = — ^r 1 , 00"= KJ ' . (74) 

dj dj 2 

The phases in the sums of the form (|59|) come into resonance, for m = 1 and hence for 
any higher m > 1, when they change by 2tt per quantum number, which happens at 
time t T 

noj t r = 2tt, t r = - — . (75) 
Too' 
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and its integer multiples*. Let time t be close to kt T , fceZ, and write t = kt T + t' where 
f < 4, so that 

fto/t' < 2tt. (76) 

Now we can estimate the general time dependent term (|59|h where g(j) is again a 
smooth classical ^-independent function of j representing a suitable combination of 
Fourier coefficients v m (j), by: (i) shifting the time variable to t f , (ii) incorporating the 
resonance condition (J75J), and (iii) approximating the resulting sum by an integral, due 
to smallness of t' ()76|) . which is a simple Gaussian: 

^#(^e im " (ftri)t A/,(M 

n 

S ^^(^)e im ("* +(nn - J * ) "' + ^ ?in -^) 2w '')(^+*') jDv) (^ n ) 

n 



= e 



V 2A ; ^ 4A ! + (W<) 2 1 

where x — ^kj* /Ti. From this calculation we deduce quantitative condition for the 
appearance and the shape of the resonance. Let Aj = ((J — j*) 2 ) 1 ^ 2 = \f^jr~ denote 
the action width of the wavepacket. Physically, we need that the coherence of linearly 
increasing phases is not lost along the size of the wavepacket, i.e. 

( := mu"tA 2 = — < 2tt. 78 

3 2A K ' 

( is precisely the coefficient appearing in the square-root prefactor and the exponential 

of the resonance profile (J77j). Indeed we see that with increasing ( the squared modulus 

of the peak of the resonance (at t' = 0) dies out as (1 + C 2 )~ 1 ^ 2 - We note that £ increases 

with the increasing order k of the resonance, since t ~ kt v , so 

TT Uj" 

( = km- — . (79) 
^ Aw' K ' 

Therefore we may get strong and numerous resonances, i.e. small (, provided either 

the second derivative to" is small, or the initial state is squeezed such that A ^> 1. For 

example, if the second derivative vanishes everywhere, uj" = 0, then the resonances 

may appear even for extended states. This is the case for our numerical model, where 

resonances can be seen also for a random state in figure |H1 

From ()77|) we read that the temporal profile in such a fidelity resonance has a shape 

of a Gaussian of effective width 

A, = (80, 
mu'Aj 

* It is interesting to note that these resonant times correspond precisely to the condition for revivals of 
the wavepacket in the forward evolution (appart from a phase-space translation) studied in Refs.[T^|. 
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modulated with an oscillation of frequency = 00*. Hence in order to feel the effect of 
the fidelity resonance, time t has to be within A t of the center of the resonance kt r . 
In the semiclassical limit, the resonance positions scale as t r oc h" 1 , while their widths 
grow only as A t oc Tr 1 ! 2 , so they are well separated. Also, with increasing order k 
the magnitudes of the peaks of the resonances decrease as ~ fc _1//2 , while their widths 
increase as ~ k, so they will eventually, at k oc h^ 1 ^ 2 , start to overlap. This will happen 
at time ~ h~ 3 ^ 2 which is smaller than ~ t 2 provided 5 < v~ l / 2 Ti. 

The resonance described in this paragraph, which will be called a 27r-resonance, 
affects all time dependent terms of (|36|37jl . but its precise shape depends on coefficients 
v m {j*)e im6 * . However, it is important to note that the fidelity can be explicitly calculated 
close to the center of the resonance, if <C A t , where the Gaussian factor of the right- 
most expression of (|77jl can be neglected. In addition we shall neglect the coefficient £ 
in (|77|) as we are particularly interested in the case of a strong resonance (<C 1. Then 
a simple calculation gives for the first moment of S t (I36|) 

(£ t > = v m (j*y m6 * (l - e im ^) , (81) 

while the second moment can be shown to be just (£ 2 ) = (S 4 ) 2 , hence the fidelity is, 
around the center of a 27r-resonance, equal to 1 within the linear response approximation 

F(t)Sl, for |t - jfct r | « A t , andC<l- (82) 

We note that such a 'flat-top' structure of 2tt— resonance is nicely illustrated in 
a numerical example in figure El where we consider a slightly modified model with 
U = exp {-i3[a(S z /S) 2 /2 + 7 {3 z /S-j*) 3 /6]), and h (j) = «j 2 /2 + 7 (j -j*) 3 /6, such 
that uj" = 7 may not be identically vanishing. 

5.3.2. 7r- and 2it /m-resonances: We note that one may obtain a resonance condition 
for time dependent term ()77j) with fixed m for even shorter time, namely for t = t T /m. 
This is trivially the case for perturbations with many, or at least more than one Fourier 
components with \m\ > 1. However, in such cases only selected time dependent terms 
of the moments (|3fi|37j) will be affected, so the fidelity will generically not come back to 
1, even in the linear response regime (JHJ and in the strongly resonant case (<1. Such 
(incomplete) resonances at fractional times (k/m)t r will be called 27r/m-resonances. 

However, we may obtain a resonant condition at t — t r /2 even for the first Fourier 
component m = 1 of the perturbation, due to taking the square of the operator E t thus 
producing Fourier number m + ml = ±2 in the last term on the RHS of eq. (J37|) . Such 
a resonant behaviour at times (k + |)t r will be called a ir— resonance. 

So for perturbations with a single Fourier mode m = ±1, or more generally with 
only odd- numbered Fourier modes m = 21 + 1, the ir— resonance can affect only the last 
term of the second moment (|37j) while it cannot affect the first moment (jSEJ)? which is 
given by its time averaged value (pjTj) . To see this, we observe that the time dependent 
parts of form g(hn)e imujt in (E t ) and (S 2 ) are proportional to Y<nD^(hn)g{Tin)(—l) mn . 
As m is an odd number and D^(hn)g(hn) is a smooth function of n this sum averages 
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Figure 4. The two-time correlation function C(t',t") © for the quantized top with 
a = 1.1, /3 = 0, 7 = 0, S = 16, and coherent initial state. 



to zero. All this allows us to again explicitly compute the linear response fidelity (|Sjl 
close to the peak in a strongly resonant case, namely 

(E t ) £ v m (f)e im9 \ for |t - (fc + \)t T \ < A t , and C < 1, (83) 
(S?) = ( E ^(/)e ime *) 2 +( £ v m (f)e^^) 2 , (84) 

m=2Z+l m=2Z+l 

4 r2 / m>0 \ 2 

F(t) — 1 7j2~ I E I«m0'*)l cos(m^*t + /3 m ) , (85) 

^ \m=2l+l J 

where (3 m are phases of complex numbers v m (j*)e ime * . So we have learned that 
the fidelity at the peak of a 7r-resonance displays an oscillatory pattern, oscillating 
precisely around the plateau value -F p i atcau with an amplitude of oscillations equal 
to 1 — -Fpi a tcau so that the fidelity comes back to 1 close to the peak of the resonance. 

Again, our numerical example illustrates such an oscillatory structure of 
7i— resonance in figure El The resonances can also be nicely seen in 'short-time' figure 
and because ( = also in the 'long-time' figure El In figure |U we depict the structure of 
7i— and 2tt— resonance as reflected in the two-time correlation function C(t',t"). Note 
that the first intersection of the soliton-like-trains for t' — t" = const and t + t' = const 
happens at t r /2 and produces a n— resonance, while the second intersection at t r produces 
a 27r-resonance. 

In analogy to the emergence of a tc— resonance as a consequence of the contribution 
from the second moment of S t , even for the first Fourier mode m = 1, we shall 
eventually obtain also fractional 27r/p-resonances at times (k/p)t T in the non-linear- 
response regimes where higher moments (Sf) contribute to F(t) ~ {exp(iH t 5/h)), eq. 
(145)1 . This is illustrated numerically in figured showing the case of strong perturbation 
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5S = 3.2 so that higher orders are important. One indeed obtains fractional resonances, 
some of which have been marked on the figure. 

5-4- Illustration in terms of echoed Wigner functions 

All the phenomena described theoretically in the preceeding subsections can be nicely 
illustrated in terms of the echoed Wigner function — the Wigner function WM s {t)\i>) (q, p) 
of the echo-dynamics. Namely, according to formula (jfi7|) the fidelity F(t) is given 
simply by the overlap of the echoed Wigner function and the Wigner function of the 
initial state. Therefore, the phase-space chart of the echoed Wigner function contains 
the most detailed information on echo-dynamics and illustrates the essential differences 
between different regimes of fidelity decay. This is shown in figure El (see [TH] ) for 
the quantized top where the Wigner function on a sphere is computed according to 
Ref.|16|. In the initial classical regime, t < t\, the echoed Wigner function has not yet 
developed negative values and is in point- wise agreement with the Liouville density of the 
classical echo-dynamics. In the plateau regime, t\ < t < t 2 , the echoed Wigner function 
decomposes into several pieces, one of which freezes at the position of the initial packet 
providing significat and constant overlap — the plateau. At very particular values of 
time, namely at the echo resonances, different pieces of the echoed Wigner function 
somehow magically recombine back into the initial state (provided ( < 1). In the 
asymptotic, ballistic regime, t > t 2 , even the frozen piece starts to drift ballistically 
away from the postion of the initial packet, thus explaining a fast Gaussian decay of 
fidelity. 

6. Semiclassical asymptotics: random initial state 

The second specific case of interest is that of a random initial state. Here we shall assume 
that our Hilbert space has a finite dimension A/", like e.g. in the case of the kicked top 
or a general quantum map with a finite classical phase space, or it is determined by 
some large classically invariant region of phase space, e.g. we may consider all states 
|n) between two energy surfaces Ei < h (hn) < E 2 of an autonomous system. In any 
case we have the scaling 

M - £ (86) 

where V is the classical d— volume of the populated action space region of interest. The 
notion of a random state refers to an ensemble average over the full Hilbert space of 
interest. So we treat the complex coefficients ip n as pairs of components of a vector on 
a 2A/"-dimensional unit sphere. In the asymptotic regime of large M , these can in turn 
be replaced by independent complex Gaussian variables with the variance 
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Figure 5. (Movies online ^H]) Snapshots of the Wigner function of echo-dynamics 
for quantized top a = 1.1, j3 = 0, 7 = with S — 200 and for S = 1.6 • 10 -3 (same as 
in figures ^JlEJi) . The upper phase-hemi-sphere is shown with j = costf € [0, 1] on the 
vertical axis and 9 = ip € [0, 2n] on the horizontal axis. From top to bottom we show: 
the initial state at t = 0, the state at t = 14 w t\ when we are around the regular 
Ehrenfest time, at t — 300 in the middle of the plateau, and at t = 100000 in the 
ballistic regime. The bar below shows the color-code of the Wigner function values. 



where ((•)) denotes an ensemble average over random states. When we write such an 
average for an operator, we actually mean 



((A)) ■.= ((mm) = 




Note that a trivial application of the pair-contraction rule (Wick theorem) yields that 
averaged fidelity is asymptotically the same as the averaged fidelity amplitude squared 

ant))) = «i/(t)i 2 )) = i«/(t)»i 2 + i = i«/(*)»i 2 . (89) 
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This means that the fidelity amplitude is self-averaging, i.e. its variance with respect 
to random state averaging is semiclassically vanishing. The same property holds for the 
fidelity itself, ((F 2 )) - ((F)) 2 = 0(1/ AT). 

6.1. The plateau: Linear response and beyond 



1 

0.98 



5S=0.32 




S=1600 



S=200 



10000 



Figure 6. Short time fidelity for a quantized top with a = 1.1,(3 — 1.4 and random 
initial state. To reduce statistical fluctuations, averaging over 20 realizations of initial 
random states is performed for S = 1600, and over 100 initial states for S = 200. The 
horizontal chain line is the semiclassical theory (|78|l . Resonances are here present due 
to the special property ui"{j) = and will be absent for a more generic unperturbed 
system. The main figure shows the case of weak perturbation SS = 0.32, whereas the 
inset shows the case of strong perturbation 5S = 3.2. 



In computing the ensemble average of the linear response formula (jHJ) we have 



to compute ensemble averages of the expressions (|36f 37J) . This is a straightforward 
application of (|%7jl for (pTTjl . and the pair-contraction rule for ({Hfijl 

r 2 sin 2 (±m ■ uj(hn)t) 

=t?E E 7^ 

■ /v n m^o sin [~m-w(nn 



T 



2 , , sin 2 ( \m ■ uf(j)t) 

TT / E KOOI 2 -^ jU, 00) 

v J m ^o sm 2 (±m ■ w(j)J 



(((S t ) 2 )) = ^((S 2 )). (91) 

We note that the expression (|9T3|) is just the classical phase space average (cx 2 ) cl where 
cr t (j, 6) = Y^=q v {3i & + ^(jW) is the classical limit of S t . We see that for a random 
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Figure 7. Short time fidelity for a = 1.1,0 = 0, S = 1600, SS = 0.32 and random 
initial state. The chain line shows the theoretical value of the plateau as computed 
from the formula (|98|l . 



state the contribution of the square of the expectation value fl9~T]) is semiclassically small 
so the plateau in fidelity is within the linear response approximation determined by the 
second moment ([90)1 . 

On one hand, let us assume that 

m ■ ^ (mod 2vr) (92) 

for all contributing Fourier components m and for all j from the classical action space 
of interest. In such a case the sin 2 in the denominator of (|9T)|) is never vanishing, so 
the integral has no problem with singularities, while the sin 2 in the numerator can be 
averaged, either over a short period of time, or over small regions of phase space for 
sufficiently long time t > t\, in both cases yielding a trivial factor of 1/2. In any case, 
the convergence timescale t 1 here is classical, and of course, does not depend on the 
strength of the perturbation t\ ~ h°5°. Hence in such a non-singular case, the linear 
response plateau of the fidelity reads 

((F(t))) = I" |U„ + 0(5% (93) 

Vi d ^ 2sin 2 (im-^))- [M) 

Going beyond the linear response approximation we again use the general formula (|43|) . 
and compute the expectation value of an operator in the random state in terms of the 
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classical phase-space average 
that is 



(2vr) d V 



d d jj d d Og(j,0), 



plateau / 



1 

V 



d d 



exp 



'z5 
ft 



(95) 



(96) 



(27r) d 

If, on the other hand, the condition ()92|) does not hold, i.e. if there exist values of 
the action j, and m e 1> d , k G Z, such that m ■ = 2nk, then the denominator of 
(j9*Uj) becomes singular and the corresponding term (in the appropriate limit) grows in 
time. However this growth stops, at least on a timescale ~ h~ x , due to a discrete nature 
of the quantum action space. At this time the value of the correlation integral, the 
plateau, will be typically so small that it cannot be described within the linear response 
approximation, so we have to employ the general formula (j45j) . The only modification 
of the general formula, with respect to a non-singular case ([96)1 . is the observation 
that the quantum phase space is in fact discrete in action so one should semiclassically 
approximate the expectation value with the sum, instead of an integral, 

1 



-y 



d d 6g(hn, 0). 



(97) 



Furthermore, the diverging terms v m / sin (m • w/2) in v m (J27|) come from the 
semiclassical approximation for S t whose quantum counterpart has matrix elements 
proportional to V mn / sin [(0 m — <f> n )/2\. As we consider perturbations with a zero 
average, V mm = 0, the diverging terms are absent in the quantum operator E t . Therefore 
to remedy the random state formula for / p i ateau (J9TI)) we simply replace an integral with 
the summation excluding the divergent terms, so the general final result reads 



((/"plateau)) 



Af 



E 



d d 

(2rt) d 



exp — v(hn, 0) 



(98) 



Again we find an excellent confirmation of our theoretical predictions in the 
numerical experiment. In the first calculations we choose the shift ft — 1.4 so that we 
have no singular frequency throughout the action space. In figure El we demonstrate the 
plateau, which in the case of random states starts earlier than for coherent states, namely 
at t\ oc h°5°. The value of the plateau can be calculated by numerically evaluating the 
integrals occurring in eq. (|94jl . for the linear response approximation, or eq. (f9T)|) in 
general. The integral over the angle 9 in the formula for the plateau ()96|) again gives a 
Bessel function so that we end up with a numerical integration over j 



((-^plateau)) 



djJ 2 5S : 



(99) 



2sin{a(j-/3)/2} / 

Observe that the random state plateau is just a square of the action-space average of 
a coherent state expression ()68|). Horizontal chain lines in figure |H1 correspond to this 
theoretical values and agree with the numerics, both for weak perturbation 5S = 0.32 
and strong perturbation 5S = 3.2 (inset). The plateau lasts up to t 2 which is for random 
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states fr- independent, t 2 ~ 1/6. Small resonances visible in the figures are due to the 
fact that the Hamiltonian is a quadratic function of the action and therefore us" = 0, so 
that the resonance condition (|75jl is satisfied also for extended states l)78j) . For a more 
generic Hamiltonian these narrow resonant spikes will be absent. In figure. 0we also 
demonstrate the plateau in the fidelity for the zero-shift case (3 = with a singular- 
frequency, = 0) = 0, where we again find an excellent agreement with the theoretical 
prediction (|98j). In this case the theoretical value has been obtained by replacing an 
integral in (jHHJ) with a sum over n (replacing j = hn) and summing over all quantum 
numbers except n = 0. Observe that the value of the plateau is much lower than in the 
case of a non-zero shift /3 = 1.4 in figure EH 

6.2. Asymptotic regime of long times 



LL 




Figure 8. Long time fidelity for random states for quantized top with a = 1.1, (3 = 0, 
for S = 200 (a), and S = 1600 (b). Here SS = 0.064 and averaging for S = 200, 
and S = 1600, is performed over 1000, and 20, initial random states, respectively. 
Heavy chain line is the theoretical asymptotic decay (|100|l with analytically computed 
prefactor (no fitting parameters) . The inset in the bottom figure shows the diffractive 
quotient between the numerical fidelity and the asymptotic formula (|101(l (chain line 
in the main figure). 

Again, after sufficiently long time t 2 ~ the second term in BCH expansion (|17|) 
will be dominating and we shall use the ensemble average of the ASI formula (J3T|) where 
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((fit))) 



y J d d j exp M 



f 

2h 



v(j)t 



(100) 



The semiclassical computation of such an integral is an elementary application of 
stationary phase method in d- dimensions, following [U] for the analogous case of a generic 
observable. The condition for the validity of stationary phase method is that tS 2 /h > 1 
which will turn out to be consistent with the assumption t > t 2 . Let j v ,r) = 1, . . . ,p 
be the p points where the phase of the exponential on RHS of (jl00|) is stationary, 
dv(jri)/dj = 0. This yields a simple result 



«/(*)» 



(2k) 



3d/2 



V 



Tl 



tr5 2 



d/2 



E 

V =l 



exp {if (jri)tr5 2 / (2h) + ip n } 



detW„|V2 



where 



jk 



d 2 v(3y. 
djkdji 



(101) 



(102) 



is a matrix of second derivatives at the stationary point rj, and v v = 7r(m + — m_)/4 
where m± are numbers of positive /negative eigenvalues of the matrix W v . Here we 
should remember that the asymptotic formula (jlUUJl has been obtained as a stationary 
phase approximation of an integral in the limit of an infinite action space. If we have 
a finite region of the action space, the stationary phase approximation of (jlOOj) gives 
additional oscillating prefactor, whose amplitude dies out as {Ti/t) 1 ^ 2 for h ^ and/or 
t — > oo, and which can be interpreted as a diffraction. This oscillating prefactor can be 
seen in numerical data for fidelity in the inset of figure |H1 

So we have found that, apart from possible oscillation due to phase differences if 
p > 1, the fidelity will for a random state asymptotically decrease with a power law 



t 



t r 



const 



h_ 



(103) 



Note that for a random initial state, the actual transition time scale t 2 , as predicted by 
eq.(|50)l. indeed does not depend on Ti. 

The above theory is again quantitatively confirmed in figure El for the quantized 
top with a = 1.1, (3 = 0, where a (single) stationary point needed for the formula 1)101)) 
had to be calculated numerically. 



7. Discussion 



In the present paper we have elaborated a semiclassical theory of quantum fidelity 
decay for systems with an integrable classical counterparts, perturbed by observables of 
vanishing time average. Such perturbations may not be generic, but provide important 
special class of perturbations which are often enforced by symmetries. 

We have found that quantum fidelity will generally, after initial decay on a short 
perturbation independent timescale t%, exhibit a saturation around a constant value 
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- the plateau, and stay there up to time t2, such that the time span of the plateau 
£2/^1 ~ l/<5 can be made arbitrary long for small perturbation 5. After the plateau, 
t > t2, the fidelity will decay as a Gaussian for a coherent initial state, or as a power 
law t~ d for random initial states, just to name the two most important specific cases, 
where the timescale of decay is generally proportional to 5~ 2 . This must be contrasted 
with the decay in the regular case of a non-zero time-average perturbation [S| where the 
decay time scales with the perturbation strength as <5 _1 . In the case of localized initial 
wavepackets in one dimension, we observe and explain the effect of echo resonances, i. e. 
the sudden revivals of fidelity at perturbation independent and equally spaced instants 
of time. 

The freezing of fidelity is a distinct quantum phenomenon, as the corresponding 
classical fidelity for the initial Gaussian wavepacket displays a power-law decay t~ d jHj 
after the point t\ where the quantum plateau starts. The classical fidelity decays on 
a time scale 5~ l no matter what the average value of the perturbation is, while the 
time scale of quantum fidelity decay drastically changes from to 5~ 2 , having V = 0. 
This increased stability of regular quantum systems to perturbations with a zero time 
average could be potentially useful in constructing quantum devices jTJj. This is even 
more so because the plateau also exists for random initial states which are expected to 
be more relevant for efficient quantum information processing. 
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